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Abstract 

The scalar form factors for the pions and kaons are calculated in SU(3) Chiral 
Perturbation Theory at order p 6 , in the isospin limit m u = = m. We find in 
general sizable corrections of 0(p & ). We use our results to obtain information on the 
1/N C suppressed low energy constants L\ and Lg as well as on two 0(p 6 ) low energy 
constants. We present some numerical results for masses and decay constants as well. 

We derive a relation for the scalar form factors analogous to Sirlin's relation for 
vector form factors. 



1 Introduction 



Chiral Perturbation Theory (ChPT) in its effective Lagrangian form was introduced by 
Weinberg pQ and developed into its standard form, for both two and three light flavours 
(referred to as SU{2) and SU{3)) formalisms, up to order p 4 , by Gasser and Leutwyler 
|21 01- The underlying assumption of ChPT is that the chiral SU{n) L x SU(n)n symmetry 
present in the QCD Lagrangian is spontaneously broken to its vector subgroup, SU{ti)l x 
SU(n)n — > SU{n)v- This breaking is due to the finite quark condensates 1 

(0\qq\0) ^ 0, for q=u,d and s. (1.1) 

The quark and gluon degrees of freedom are integrated out and replaced by pseudo-scalar 
fields representing the lightest pseudo-scalar mesons triplet (577(2)) or octet (S77(3)), 
these are the Goldstone bosons resulting from the spontaneous chiral symmetry breaking. 
Introducing non-zero quark masses renders the chiral symmetry only approximate. This 
results in masses for the mesons. 

ChPT is thus an effective theory of QCD, perturbatively expanded in powers of the 
relevant momentum and the quark masses, and is non-renormalizable. Loop corrections 
produce ultra-violet divergences which require the introduction of new counterterms and 
associated low energy constants (LEC) at each order in the parameters expansion. The 
number of new operators is, however, finite at each order, allowing ChPT to be predictive. 
Recent introductions to ChPT can be found in Ref. j3]. 

A basic set of calculations in S77(3) ChPT at order p 4 was performed in [13 03 and used 
to determine the LEC at order p 4 , the L\, to that order. L\ and U % were assumed to be 
negligible there using the large N c limit of QCD This assumption was studied at order 
p 4 in [7] by considering scalar form factors and a scalar two-point function. 

The state of ChPT in the meson sector now is such that many calculations have been 
pushed to the next order, p 6 or next-to-next-to-leading (NNLO) or two-loop order. There 
are two main reasons to explore the two-loop region of S77(3) Chiral Theory. First, the 
accuracy of experimental results has improved so this accuracy is now routinely needed 
and we need to determine all LECs to the required precision. Second, spontaneous chiral 
symmetry breaking of QCD does not have to be driven by the quark condensates, Eq. 
(jl.ljl . An alternative bookkeeping method in the Lagrangian expansion might be neces- 
sary, leading to a Generalized ChPT (GChPT), see Ref. [S] and references therein. The 
combination of the two-loop calculation of ^-scattering in S77(2) ChPT [31101, together 
with the dispersive Roy equation analysis [TT], lead to sharp predictions for the tttt scat- 
tering lengths ^21 confirmed by the E865 experiment leading to the conclusion that 
SU{2) ChPT is indeed driven by the quark condensate [Hj. See also [HI [13112] for relevant 
discussions. In the case of SU(3) ChPT the possibility of a small quark condensate is not 
ruled out 0[T71[TB|. The pion scalar form factor has been evaluated to order p % in S77(2) 
ChPT in ^1 after earlier dispersive work [20J. 

1 See below for the alternative view of generalized ChPT. 
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The two-loop calculations in SU(3) ChPT of the masses [2B 122] indicated the possi- 
bility of sizable NNLO corrections. This behaviour persisted when the calculations also 
included Ku j22l Elj such that the L\ could be determined with p 6 accuracy under similar 
assumptions as the p 4 determination In j21] a search was done for possible values of 
L\ and L r 6 that would lead to acceptable convergence of all quantities considered there. In 
[2*5*] effects of isospin violation in the masses were taken into account also to order p 6 as 
well as well as the new E865 data of J3] ■ 

This paper consists of two major parts. In the first part we calculate the scalar form 
factors of pions and kaons to order p 6 in SU (3) ChPT and show some numerical results 
for the sizes of the various contributions to the different form factors. We therefore first 
present the framework of SU(3) ChPT in the isospin limit and including external scalar 
fields in Section |2] and sketch the calculation of the scalar form factors in Section |3] The 
analytical results of this calculation are given in Section 14.11 and App. El In Section 14.31 
we present a series of plots allowing the reader to judge the importance of the different 
parts of the calculations. One result of more general interest is that the curvature of the 
scalar form factor in decays [2E] can in principle be determined from the quantities 
considered here. Unfortunately as discussed at the end of Sect. I5.5l the numerical accuracy 
obtainable is rather low. 

Those results allow us then to proceed to the second major part of this paper. A study 
of the form factors allowing to relax some of the assumptions made in the previous L\ 
determinations and the effects of these assumptions on the scalar form factors. We made 
use of dispersive analysis, Section 15.11 to extract experimental information. This allows us 
to obtain extra numerical results for some of the LEC's of the theory, Sections 15 . 31 and 15.51 
We also present some updated results on masses and decay constants in Sect. 15.61 

A third result is a relation between the scalar form factors that has only second order 
corrections in quark mass differences. This relation, Eq. (|3.18|) is proven in general in 
App. \M and is valid for all t where an expansion in mass differences is valid. It is the 
analog of the relation derived by Sirlin for the vector form factors [2Z] • 

2 The ChPT Lagrangian 

This section contains a very short overview of ChPT and serves also to define our con- 
ventions. More detailed introductions can be found in 4|. The Lagrangians of ChPT are 
ordered in the p counting, powers of momenta (0(p)), quark masses or external scalar or 
pseudo scalar fields {0{p 2 )) and external vector or axial- vector fields (0(p)). We keep to 
the three flavour case here. 

The Lagrangian for the strong and semi leptonic mesonic sector to NNLO can be written 

as 

£ = £ 2 + £ 4 + £ 6 , (2.2) 
where the subscript refers to the chiral order. The lowest order Lagrangian 

£ 2 = S(«X + X+>- (2-3) 
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The mesonic fields enter via 

/ %M 
u = exp — 

\F0V2 



M 



\ 



7i 

K~ 



7T+ K+ \ 



(2.4) 



and the quantity introduces the external vector (t> M ) and axial-vector (a M ) currents 



u, 



iiu^d^u - d^uu ] - iu ] r^u + iul^u ] ) , = ty, - (+)< 



while the scalar (s) and pseudo scalar (p) currents are contained in 

X± = ± u^it , X = 2S (s + ip) . 

The p 4 or NLO Lagrangian £ 4 was introduced in Ref. [3] and reads 

+L 3 (/n/)i iy ) + L 4 (m%)(x+) 
+L 7 (x-) 2 + L 8 (x+X-) 

-iL 9 {F^uu^u u u ] + F^u"u v u) + L w (F*UF L ^U j ) 
+H 1 (F^ R + F^ L ) + ^(x 2 , - x 2 >/4 . 

The L 9 and L w terms introduce also the field strength tensor 

Ku R) = Wr)u - dj(r), - 1 [l(r)„ J(r)„] . 



(2.5) 



(2.6) 



(2.7) 



(2.8) 



The two terms proportional to H\ and H2 are contact terms and do not enter physical 
amplitudes. In the present case, we keep only the relevant scalar current interactions from 
which we extract and separate the quark masses contribution: 



+ 5, p = l. 



0. 



(2.9) 



We quote the schematic form of the NNLO Lagrangian in the SU(3) case 

^6 = Cj Oi 



(2.10) 



i=l,9 



and refer to [28 , where this was first constructed after the earlier attempt of [29], for the 
full expressions. The last four terms are contact terms |28j . 

All ultra-violet divergences produced by loop diagrams of order p 4 and p 6 cancel in the 
process of renormalization with the divergences extracted from the low energy constants 
L^s and Cj's. We use here dimensional regularization and the standard modified minimal 
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subtraction (MS) version used in ChPT. An exhaustive description of the regularization 
and renormalization procedure including a description of the freedom involved can be found 
in Ref. [E] and [H]. 

The subtraction of divergences is done explicitly by 



and 



a 



2(<2-4) 



F 2 



1 ( 2 ) A 2 



where c and A are defined by 



(2.11) 
(2.12) 



lnc= --[ln47r + r'(l) + l] 



(2.13) 



A 



16vr 2 (rf-4) 



The coefficients T,-, T'P and T ( /' are constants while the T) 1 " 1 's are linear combinations 



,(2) 



(2.14) 



of the Ll(n) 's. Their explicit expressions can all be found in jSUj where they have been 
calculated in general. 

Scalar form factors can be calculated using functional derivatives w.r.t. the external 
scalar field s |3j- 



3 The scalar form factors: definitions and overview 
of the calculation 

3.1 Definitions 

The scalar form factors for the pions and kaons are defined as follows: 

(M 2 (p)\q t q 3 \M 1 (q))=F^(t) (3.15) 

with t = (p — q) 2 and i, j = u,d, s being indices in the flavour basis and Mi, M 2 a meson 
state with the indicated momentum. 

In the present case of isospin symmetry m u = = rh, the various pion scalar form 
factors obey 

F§(t) = ^\t) = 2i^ V °(t) = 2i^ + (t) = 2F^ + (t) = 
= -2v / 2F^ + (t) = 2v / 2F M ^-(t), 
Fl = FZ + * + = ■ (3-16) 
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The kaon currents are related by the following rotations in flavour space 

Flit) es F« +K+ V = F«°«\t), 

Fs K d (t) = Fr^ + (t)=Fr x °(t), 
- F s r K+ (t)=F s K s ° K °(t) 1 

Ft H) = F S K ^ it) = V2F^° it) , 

F£ q {t) = + (3-17) 

The other scalar form factors can be related to these using charge conjugation and time 
reversal. In ()3.16|) and (j3.17j) we have also given the notation we shall use for the form 
factors in the remainder. 

The scalar form factor Fg(t) is proportional to the form factor fo{t) used in 
decays. 

The scalar form factors can be shown to obey a relation similar to the Sirlin [27] relation 
for the vector form factor 

Fg(t) - 2Fl(t) - 2F«(t) + 2Ff,{t) - 4Ff *(t) = O (K - ™?) • (3-18) 

The proof of this relation is in App. 1X1 

The values at zero momentum transfer are related to the derivatives of the masses w.r.t. 
to quark masses because of the Feynman-Hellman theorem (see e.g. |Hj) 



^(0) = ^, FU0) = J-m 2 



dm 7r ' * s w dm, 
d n d 



7T ) 



^(0) = ^. (3.19) 

This also shows that we can expect large corrections at t = 0. The argument is fairly 
simple. If 

m\ « B m s + (3iB m s ) 2 + 7(5 m s ) 3 , (3.20) 

then 

F£(0) « B + 2(3B m s + 37(5 m s ) 2 . (3.21) 

So we see that in the scalar form factors the relative p 6 corrections can get enhanced by 
factors of order 3 compared to the masses. 



3.2 Diagrams 

The relevant diagrams for the present case are identical to those involved in the electro- 
magnetic form factors [31] with the electromagnetic current replaced by the scalar one. We 
list the diagrams appearing at each order but refer to j2H EH] for a deeper discussion of 
the checks to perform and the renormalization. 
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a) b) c) d) 



Figure 1: Order p 2 and p 4 . 

3.2.1 Leading and next-to-leading orders 

The diagrams are shown in Fig^ The black square represents an external scalar insertion 
from £ 2 - diagram a) - and the crossed square a scalar insertion form £ 4 - diagram d). 
Also included, are the tadpole contributions with two possible insertions. 
The lowest order contributions are 

F§(t) = 2B , Fl(t) = 0, 
F«(t) = B , F&{t)=B 

Ft(t) = B . (3.22) 

3.2.2 Next-to-next-to-leading order 

There are four distinct topologies involved in order p 6 diagrams. Firstly, the two loop 
diagrams, built exclusively on 0(p 2 ) vertices, are represented in Fig El There again does a 
black square represent a scalar interaction. 

Secondly, one loop diagrams can contribute, with one vertex built on £ 4 , possibly 
including a scalar insertion. These are represented in FigEl 

There are also sunset integrals and irreducible contributions to include through the 
non-factorizable diagrams of Fig |3J See [21] and jST] for a treatment of sunset integrals 
and irreducible two loop integrals. 

Of course, one must finally include the £ 6 tree contribution, Fig|5] 

4 Analytical and First Numerical Results 

4.1 Analytical expressions - 0(p 6 ) LEC's 

We present here the dependence of all the scalar form factors on the LECs of order p 6 . It 
can be easily checked that the relation (j3.18|) is satisfied. 

(7r|uii + dd\n) = ^{ml m\ (-16C[ 3 - 8C[ 5 + 8C[ 6 + 8C r 20 + 72C r 21 + 16C r 32 ) 
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c) d) e) f) 



Figure 2: The two-loop corrections diagrams. 




c) d) e) 



Figure 3: One-loop diagrams involving an 0(p 4 ) vertex. 




a) b) 



Figure 4: Irreducible two-loop diagrams. 
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Figure 5: The Cq contribution. 

+ ml t (4C[ 2 + 8C[ 3 + 2C r u + 3C[ 5 + 4C[ 6 + 2C{ 7 + 2C r u + 2C r 36 ) 
+ ml (-12C[ 2 - 16C[ 3 - 6C[ 4 - 8C[ 5 - 14C[ 6 - 6C[ 7 + 18C[ 9 
+ 30C 2 r + 30C^ + 12CJi + 16C 3 r 2 ) + 2m\ t C[ 5 

+ m% (SC r 16 + 8C r 20 + 2AC r 21 ) - t 2 (C[ 2 + 2C{ 3 ) } . (4.23) 
H^tt) = ^{mlml (-16C[ 6 + 16^ + 48^) + m*t (8C[ 3 + 2C[ 5 - 4C[ 6 + 2CJ.) 

7T 

+ ml (-8C[ 3 - 4C[ 5 + 8C[ 6 + 24C^ + 8C&) + 8m 2 K t C r 16 - 2 1 2 C[ 3 (}.24) 
(K + \uu + dd\K + ) = ^I + m l m 2 K (-8C r 13 -AC r 15 -8C r 16 + 16C r 20 + 72C r 21 + 8C r 32 ) 

7V 

+ m 2 t(2C[ 4 + C[ 5 + 8C[ 6 -2C[ 7 ) 
+ ml {-2C r u - 6C[ 6 + 2C[ 7 + 6C[ 9 + QC r 20 + 6C r 21 ) 
+ m\ t (4C[ 2 + 16C[ 3 + 6C[ 5 + 4C£ 7 + + 4C£,) 
+ (-12C[ 2 - 32C[ 3 - 4C[ 4 - 16C[ 5 - 8C[ 6 - 8C[ 7 + 12C[ 9 

+ 40C£ + 120^ + 12C r 31 + 32C r 32 ) - t 2 (C[ 2 + 4C[ 3 ) } (4.25) 
^L{+mlm 2 K {-8C{ 3 + 8C[ 4 - 4C[ 5 + WC[ 6 - 8C{ 7 - 24C[ 9 - 8C r 20 

IT 

+ 48C^ + 8C r 32 ) + ml t (-2C[ 4 + C{ 5 - 4C[ 6 + 2C[ 7 ) 
+ ml {-2C r u - 6C[ 6 + 2C[ 7 + 6C[ 9 + 6C 20 + 6(7^) 
+ m£ f (4C[ 2 + 8C[ 3 + 4C[ 4 + 4C[ 5 + 8CJ. + + 2CJ,) 
+ (-12C[ 2 - 24C[ 3 - 12C[ 4 - 12C[ 5 - 24C[ 6 + 36C[ 9 + 48C 2 r 

+ 72C r 21 + l2C r 31 + 2AC r 32 ) - t 2 (C[ 2 + 2C[ 3 ) } . (4.26) 
^{ + ml m 2 K (-2C[ 2 - 12C[ 3 - 6C[ 5 + 8C[ 6 - 4C[ 7 + 4CJ, + 2AC r 21 

T 

+ 4C& + 12C; 2 + 2C r u ) + ml t {2C[ 2 + C[ 5 + 2C{ 7 + C r u ) 



(K + \ss\K + ) = 



(ti~\su\K ) = 
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+ m 4 (-5C[ 2 - 4C[ 3 - 2C[ 4 - 2C[ 5 - 6C[ 6 - 2C[ 7 + 6C[ 9 + 10C r 20 

+ 6C^ + 4C3j + 4C£ 2 - CJJ + m 2 , t (2C[ 2 + 2C[ 4 + 2C[ 5 + C£J 

+ m 4 K {-hC{ 2 - 8C[ 3 - 4C[ 4 - 4C[ 5 - 8C[ 6 + 12C[ g + 16C£ + 24^ 

+ 4C 31 + 8CJ, - C 34 ) - t 2 C[ 2 } (4.27) 

4.2 Loop corrections 

We present only the results for Fg. The expressions for the others are obtainable on request 
from the authors. Form factors including external etas can be calculated as well if required. 
The order p 4 results are in agreement with those of [Sj. The form given in App. El is the 
one that corresponds to the p 6 expressions we use and which can also be found in App. [Bl 

4.3 First Numerical Results 

In this part we present some plots of the sizes of the various corrections for the case with 
C\ = 0, we have used the L\ values of [2*3] . fit 10, and the neutral pion, neutral kaon and 
physical eta masses as well as F n = 92.4 MeV and /i = 770 MeV. These plots are included 
here to show the relative sizes of the many contributions to the scalar form factors we have 
calculated. These corrections are in many cases large as expected from the argument given 
at the end of Section Kill 

In Fig. |HK we have plotted the loop contributions to the pion scalar form factor nor- 
malized to 2B . The order p 2 is then one. The two-loop contribution in this case is sizable 
but not enormously so. The loop contributions to the pion strange scalar form factor, 
Fg s (t)/B , are plotted in Fig. Eb- The corrections here are moderate. Note that the imagi- 
nary part vanishes at order p 4 and remains small at p 6 . This was expected since it requires 
a kaon or eta loop when the LECs are set to zero. 

In Fig. [?li we have plotted the loop contributions to the kaon light quark scalar form 
factor normalized to Bq. There are large two-loop contributions in this case for the real 
part while the imaginary part has only modest corrections. This agrees with the imaginary 
part calculated previously in [32] • 

The kaon strange scalar form factor, F^it) / Bq, is plotted in Fig. Efc>- The corrections 
here are also large compared to the order p 2 result which is 1. The imaginary parts vanish 
at order p 4 and remain very small at order p 6 . 

We will do a more extensive analysis in the next section but here we want to show the 
effects of the diagrams involving L\ separately. We use here the L\ of fit 10 of (23, the 
Cl(fj, = 770 MeV) = and the other input parameters as given above. We also present 
the dependence on L\ by changing it from zero to L 4 (/i = 770 MeV) = —0.0003 keeping 
the others at the values given by fit 10. 

In Fig. |HK we plotted the effect of the L\ on Fg (t) / (2B ) compared to the case with the 
LI = 0. The curves are for the real parts at order p 4 and order p 6 . The curves are labelled 
respectively Fit 10, L r { = and L r A ^ 0. Notice that the effect of such a small U A is fairly 
large. 
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loops only L- = C| r = 



loops only L- = C- = 




-0.2-0.15-0.1 -0.05 0.05 0.1 0.15 0.2 0.25 

t [GeV 2 ] 

fa) 



CO 




-0.2-0.15-0.1 -0.05 0.05 0.1 0.15 0.2 0.25 

/2i 



t [GeV £ ] 

(b) 



Figure 6: The various contributions to (a) Fg(t)/(2B ) and (b) Fg a (t)/B as a function of 
t for the case with L\ = C\ = 0. Plotted are the real and imaginary parts of the p 4 and 
the p 6 contributions separately. 




Figure 7: The various contributions to (a) Fg q (t)/B and (b) Fg{t)/B Q as a function of t 
for the case with L\ = CI = 0. Plotted are the real and imaginary parts of the p 4 and the 
p 6 contributions separately. 
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U Fit 10 C[ = 



Li Fit 10 C-=0 
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m 

in 
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t [GeV 2 ] 

(a) 



0.15 
0.1 
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-0.05 
-0.1 
-0.15 



I I 


Rep 4 L i= 




Re p 4 L.4*0 - 




Rep 6 Fit 10 




•VX Rep 6 L i= _ 




- . Re p 6 L 4 *0 




% ' - 


1 1 





-0.2-0.15-0.1 -0.05 0.05 0.1 0.15 0.2 0.25 
,2-, 



t [GeV ] 

(b) 



Figure 8: The effects of the L\ on the various contributions to (a) Fg(t)/(2B ) and (b) 
Fg s (t)/B as a function of t for the case with CJ = 0. The curves are for the real parts at 
order p 4 and order p 6 . The curves labelled respectively Fit 10, L\ ^ and L\ = are for 



the standard values of the U, 



the same values but L r A = —0.003 and with all LJ = 0. 



The pion strange scalar form factor, Fg s (t)/Bo is shown in Fig. |Hb- At order p A only 
L r A contributes but it has a large effect. Lg contributes as well but is assumed zero. So 
there is no difference for this case for Fit 10 and the Lj = 0. Notice that the effect of the 
p 6 corrections is larger than the change due to this value of L\. It is therefore necessary 
to include the p 6 effects before drawing conclusions on the values of L\ and L\. 

We can present similar plots for the L[ dependence of the kaon scalar form factors but 
the conclusions are qualitatively the same as for the pion form factors. 



5 Numerical Analysis 
5.1 Dispersive Inputs 

There is no direct experimental measurement of the scalar form factors but information 
can still be extracted from elastic scattering experiments. We recall here the principal 
results of dispersive analyses and refer to [UHl and [7] for a more thorough treatment. An 
analysis within a unitarized model approach can be found in |34j . 

5.1.1 The Muskhelishvili-Omnes (MO) problem 

We consider first the one channel case, i.e. where only pion interactions are considered. 
The problem of finding the functions satisfying the pion form factor's analytical properties 
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is referred to as the Muskhelishvili-Omnes problem. Provided some convergence conditions 
on the form factor and the phase shift S^s) are satisfied, the solution is of the form: 



F w (s) = P(s) exp 



as - 



7T 



s' 0' - s) 



(5.26 



Watson's theorem then relates the phase shift to the scattering phases 5j, here 5 W = Sq. 
This can be extended to two or more channels. We restrict here to the two channel case, 
7T7T and KK in the isospin and angular momentum zero case, 1 = S-w&ve. The MO 
problem becomes a system of coupled equations of the two independent contributions 

F l {s) = (0|X |7T7r) and F 2 (s) = -^(0\X o \KK) , (5.29) 

v3 

where X Q is an isospin zero operator, e.g. (uu + dd) or ss. The interesting feature of the 
MO solution is the possibility to write the general solution using two independent sets of 
solutions {F±\ }i=i,2 fulfilling the initial conditions 

F^(0) = 5 lv (5.30) 

and this independently of the exact form of the operator X . The solution is then built on 
the value of the form factors at the origin: 

= F 1 (0)Fi 1 \s) + F 2 (0)Fi 2 \s) 
= F 1 (0)F 2 {1 \s) + F 2 (0)F 2 i2 \s) 




5.1.2 The solutions 

The method and coding of the solutions were borrowed from Moussallam ESI, where 
the solutions are obtained by solving the linear system of equations - equivalent to ()5.28|) 
for the two channel case - given by the discretized form of 



ReFi(s) =PV I ds , lmF ^ S > (5.32) 

J 4 m 2 S S 

and 

ImFi(s) = 2i ^j(s)T i j(s)F j (sy (5.33) 

3=1,2 

where 




(jj(s) = \ 1 6(s — Am-) , rrij = m^, m K (5.34) 



and Ty are the T-matrix elements of the needed tttt, KK scattering channels. Three 
T- matrix models of scattering are used, from Au et al. j^Hl, Kamihski et al. jH7j and 
Ananthanaryan et al. [S3], and we refer again to [7J El] for details and differences. The 
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F 1 (0) 


F[(0) 


1/2^(0) 


F 2 (0) 




1/2F»(0) 


Au« 


1. 


2.33 


9.94 


0. 


0.92 


5.31 


Au( 2 ) 


0. 


0.31 


1.03 


1. 


0.78 


0.98 


l m (i) 


1. 


2.43 


10.28 


0. 


0.75 


3.49 


l m (2) 


0. 


0.27 


0.93 


1. 


0.84 


1.00 


Abm« 


1. 


2.45 


10.61 


0. 


1.03 


5.36 


Abm( 2 ) 


0. 


0.21 


0.72 


1. 


0.81 


0.98 



Table 1: Results of the numerical computation of the Omnes solutions for the different 
models. The superscript refers to the respective inequivalent sets of solutions. The first 
two solutions are from Ref. jZj, the last from |35j . 

results of our calculation for the different derivatives of the canonical solutions using the 
programs of EH] at the origin are presented in Table We have only used the last 
solution, since they use a better physical model for the tttt to KK amplitude and newer 
7T7T phases. 

Similar work exists for the Kir form factor, Fg n (t), in [HE]. We have not obtained their 
results in a form we can use in the same way. 

5.2 Other experimental inputs and Values of the L[, i ^ 4, 6 

The other experimental inputs are the same as those used in [2H] corresponding to fit 10 
there. This differed from the main fits reported there by using the then preliminary data 
of H3]. 

We would like to study the dependence on L\ and U % as well. A study was already 
performed in [21] where we varied the input assumptions of L\ and U % and then refitted 
the other L\. We have redone this now with the data of 13. included, so we use all the 
inputs described above with as input values for L\{m^) and Lg(m p ) values on the grid 

\{fL\ = -4,. ..,6 and 10 4 L r 6 = -3, . . . , 6 . (5.35) 

In all these cases we have used the resonance saturation contributions to the constants 
of order p 6 as derived in [21] and those for the decay constants and the masses as derived 
in [21] but we have chosen the value of d m as defined in [21] to be zero. This is also the 
choice made in [21]. The reason for this is that the estimate of the parameter of [21] was 
very uncertain and the naive value obtained there led to p 6 contributions to the masses 
which were enormous. 

Fits of a quality similar to the one with L\ = U & = could be obtained for most points 
on the grid (|5.35|) . All those satisfying 

Lg < O.6L4 + 0.0004 . (5.36) 
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<ji|uu + 3d|it> <K|uu+9d|K> 




Figure 9: (a) The result for Fi(0) ud = Fg(0)/B as a function of L\ and L\. Also shown 
is the lowest order value 2. (b) The result for F 2 (0) u d = 2/v / 3i 71 J^(0)/i?o as a function of 
L\ and Lg . Also shown is the lowest order value 2 / \f?> . 

had a value of x 2 below 1. The precise definition of the \ 2 can be found in |2"4*1 125j. The 
data used up till now do not allow for a determination of L\ and Lg. But note that while 
these cannot be determined, a given assumption on their value leads to correlated values 
of all the other L[ in order to obtain a decent fit to the data. This effect is taken into 
account in the next subsections. 

In the numerical results quoted below we have used the neutral kaon and pion masses 
and a subtraction scale /i = 770 MeV. 

5.3 Variation with L\ and Lg of the scalar form factors at t = 0. 

It should be kept in mind that the values of the other L\ used change as well with the 
values of L\ and L\ in accordance with the fit to the other experimental values as described 
in Section I5~21 The C r -(m p ) contributing to the scalar form factors are set to zero here. 

We first show the results for Fg(0)/B in Fig. E^. We also show the lowest order 
expectation of 2 for comparison. Only the points with good fits are shown. 

We show similarly in Fig. the value of F 2 (0) u d = 2/ \/3 Fg / B together with the 
lowest order value of 2/\/3. Notice that in both cases the corrections are large for most 
values of L\ and Lg but reasonable in the region of the axis Lg ~ L\ — 0.0003. Requiring 
the corrections to be of a reasonable size thus leads to constraints on L\ and L\. 

We have also studied the ss current. The results for the pion case are shown in Fig.fTUk. 
The result for F 2 (0) s = 2 A/3 Fa JO) / Bn is shown in Fig.lTUb. Notice that the strange quark 
content in the pion remains small for the entire studied region. 
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<^|ss|tu> 



<K|ss|K> 




(a) (b) 

Figure 10: (a) The result for Fi(0) s = Fg s (0)/B as a function of L\ and U % . Also shown 
is the lowest order value 0. Notice the strange quark content of the pion remains small 
over the entire studied area, (b) The result for F 2 (0) s = 2/V3F^ s (0)/B function of 
L r 4 and L\. Also shown is the lowest order value 2/\/3. 



At this point we cannot determine the values of L r A and U 6 but there is also a tendency 
in those two cases for the corrections compared to the lowest order to be small for positive 
L\ and a correlated value for U & along U & ~ L\ — 0.0003. 

5.4 Slopes and the value of L\ and L\ 

We can use our results to get at the form factors away from zero in two different ways. 
First we can of course simply use our full ChPT calculation to calculate them to order p 6 
using the input values of L\ and L\ and the resulting fitted values of the other L\. 

Second, from Sect. 15.11 we can calculate the form factors at t ^ given the values at 
zero, this uses then all the dispersive constraints. 

We then require that both methods give the same results and obtain in this way con- 
straints on the values of L\ and F\. The main effect of this comparison can already be 
studied by looking at the radii only since for the curvature we expect a larger effect from 
the C\. We still assume here that the C[(m p ) = for those contributing to the radii. 

The results for the pion scalar radius, 

< r 2 >^= -4tt^WIw , (5-37) 

are shown in Fig.lTTk. We have normalized here to the value at zero. That leads to singular 
predictions when Fg(0) vanishes, but away from that region the dispersive prediction is 
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(a) (b) 

Figure 11: (a) The result for the scalar radius < r 2 >f as a function of L\ and Lg, with 
the dispersive estimate (abm) and the ChPT calculation, (b) The result for the kaon light 
quark scalar radius < p 2 >f- as a function of L\ and Lg, with the dispersive estimate 
(abm) and the ChPT calculation. 



quite stable and is entirely predicted by the values of the tttt phase shifts. It is fully 
dominated by the first independent solution of the MO equations. The ChPT prediction 
is quite dependent on L\ and Lg The two agree for positive values of L\ and a correlated 
value of Lg. 

The results for the kaon light quark scalar radius, 

are not as easily shown, as the calculated value of Lj^(O) runs through zero in the relevant 
region, producing strong effects in both the dispersive and ChPT results. We plot instead 
the scalar radius normalized to the lowest order result 

< P 2 >s q = j- o j/s q (t)\t=° ■ (5-39) 

This is shown in Fig. II lb . The dispersive prediction is less stable here. The relative 
strength of the two canonical MO solutions varies much more than in the pion case. This 
has a strong effect since both solutions have a very different radius, as can be seen by 
comparing the first derivatives columns in Table d 

Notice that the agreement between the chiral and dispersive radii tend to be in the 
same region of L\ and Lg where the corrections to the lowest order results for the form 
factors at t — are fairly small, as discussed in the previous subsection, taking into account 
that the results are less reliable in the kaon case. 
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<K|ss|K> 




(a) (b) 

Figure 12: (a) The result for the pion strange quark scalar radius < p 2 >f s as a function 
of L\ and Lg, with the dispersive estimate (abm) and the ChPT calculation, (b) The result 
for the kaon strange quark scalar radius < p 2 >|- s as a function of L\ and Lg, with the 
dispersive estimate (abm) and the ChPT calculation. 



The strange quark scalar current can be studied in exactly the same way. The pion 
strange scalar form factor has a very small value at zero, we normalize therefore to the 
value B and plot 

<P 2 >Ss=-^J t Fl(t)\t=o. (5.40) 

The results are shown in Fig. 112b both for the dispersive and the ChPT result. Again the 
dispersive result shows a significant variation with L\ and Lg since both solutions of the 
MO equations are important and their contributions vary significantly depending on the 
values of L\ and Lg. 

The prediction for the kaon strange quark radius 

is stable for most values of the L£, Lg except for a region near the top of the L\ and 
the bottom of the L r R range where F£(0) is near zero. The reason is that for most input 
values the second MO solution dominates. The two predictions never agree, the dispersive 
estimate is positive and the ChPT one is negative. We plot instead 

< P 2 >ss= j- Q J t Fss(t)\t=o , (5.42) 
in Fig. IT2"b . For this quantity there is a small region of overlap between the dispersive and 
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the ChPT estimates but it is where Fg s (0) is near zero and we can have a large sensitivity 
to effects from the G\. 

It is somewhat difficult to get a final conclusion about L\ and L\ since the effect of 
the p 6 constants has to be evaluated. This requires a general study of effects in the scalar 
sector which goes beyond the scope of this paper. 2 The main constraint from the pion 
scalar radius is 

Ll&L r 4 - 0.00035. (5.43) 

To be precise, this is the curve where the two surfaces shown in Fig. ITTk intersect. If 
in addition we require that the values of the scalar form factors at zero do not deviate 
too much from their lowest order values we obtain that L\ should be in the range from 
0.0003 to 0.0006. The latter requirement means staying close to the intersection of the two 
surfaces in Fig. Ek- 

We have shown several numerical results in Table El for grid points in this region. 



5.5 Curvatures and the value of C[ 2 and C[ 3 

Since the values at zero determine the form factors away from zero, we can also compare 
higher order terms in the expansion in t. Because of the large corrections to the values at 
zero coming from order p 6 , there is a rather large uncertainty in the resulting values. We 
have therefore not done a full error analysis. 
We define the pion scalar curvature via 

c -?sW^ F|(t)U - (5 - 44) 

The results both from the dispersive and the ChPT analysis are shown in Fig. 113b . It 
should be remarked that, precisely as in the two flavour case analyzed in ^U] a large part 
of the curvature is due to the loop effects and the contribution from the p 6 constants is 
not the dominant one. This can be seen by looking at the scale in Fig. IT3*k . The difference 
between the ChPT prediction with the C|" = and the dispersive result is small. 

The same analysis can be performed for the other form factors. We again run into the 
problem of normalizing the curvature, just as we had for the slopes. For the rest we use 
thus quantities normalized to Bq. We show the pion strange scalar curvature 

in Fig. [T3*b . The variation of the dispersive result is large here again due to the fact that 
both MO solutions contribute with a strong variation in relative strength. 
Analogously to (|5.45j) we define 

1 1 d 2 

ls ~ 2K2dr* Fs{t)lt = o > 



2 The analysis might be doable by extending some of the lines of work of |39| . 
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(a) (b) 

Figure 13: (a) The result for the pion scalar curvature function of U A and Lg, 

with the dispersive estimate (abm) and the ChPT calculation, (b) The result for the pion 
strange scalar curvature 7^ s as a function of L\ and Lg, with the dispersive estimate (abm) 
and the ChPT calculation. 

lid 2 

but we have not included plots of these quantities. 

The curvatures can now be used to estimate combinations of the C\. We get 

n r i or r K ( ndisp tt ChPT,CV=0\ __ / K disp KChPT,Cl=0\ 

o 12 + zo 13 — — \c s c s j — — y~f Ss Iss J > 

r r K ( ndisp 7rChPT,Cl=0\ 

cr 2+ 4cr 3 = (5.47) 

From the figures IT3*k and fTTTb it is obvious that the extracted values will depend strongly 
on the input values used for L\ and Lg, since the difference between the surfaces in those 
figures gives the values of the C\. On the other hand, we see that there are several ways 
to extract the same quantities. It turns out that these tend to give similar values in the 
region of L\ preferred in the previous subsections. Away from these regions, the obtained 
values for the C\ are much larger than naive order of magnitude estimates. 

We have collected in Tabled various results for four sets of L\. The first one is exactly 
the same as fit 10 of Ref. and all quantities coincide with the one reported there. The 
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Set 




Fit 10 


A 


B 


C 


10 3 • L\ 




0.0 
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10 3 • Lg 
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Fff(0)/B n CChPT C r = 


(Ti 


u.uzu 
1 QzL 

1.77 


n nru 

l.OU 

1.35 


n nn9 

— u.uuz 
1 1 9 

1.17 


u.uuo 
1 ^1 

1.45 


S ^iiii j 

2 n,ChPT,Cr=0 2 


) 


n ri 7 


U.O^fc ( 


n (Si n 

u.ozo 


n fii 4 

U.U14: 

u.ooo 


10 (Cl2 + 2C T 3 ) 1 C S\ 
1U (U 12 + ^13) [1s\ 

10 5 (CT 2 + 2CT 3 ) h#J 
10 5 C[ 3 [ 7 SJ 




—2.6 


—0.56 


n r r 

0.55 


—0.71 




— o.o 


— u.oo 




— U. 1 o 




-0.55 
-0.56 


0.11 
-0.02 


0.33 
0.15 


0.15 
0.03 


10 5 (C[ 2 + 4C[ 3 ) [ 7 fJ 




-4.1 


-0.27 


0.99 


-0.08 






-1.5 


-0.52 


0.26 


-0.78 


10 5 C[ 2 [c5, 7 fJ 




-1.1 


-0.84 


0.12 


-1.3 


F (MeV) 




87.7 


63.5 


70.4 


71.0 


F (4) 

7T + 




0.136 


0.230 


0.253 


0.254 


F (6) 

7T + 




-0.083 


0.226 


0.059 


0.048 






0.169 


0.157 


0.153 


0.159 






0.051 


0.063 


0.067 


0.061 


u (2) 




0.736 


1.005 


1.129 


0.936 


« (4) 




0.006 


-0.090 


-0.138 


-0.043 


U (6) 




0.258 


0.085 


0.009 


0.107 


(2) 




0.687 


0.938 


1.055 


0.874 


(4) 




0.007 


-0.100 


-0.149 


-0.057 


(6) 




0.306 


0.162 


0.094 


0.183 


(2) 




0.734 


1.001 


1.124 


0.933 


(4) 
/V 




-0.052 


-0.151 


-0.197 


-0.104 


(6) 
/V 




0.318 


0.150 


0.073 


0.171 



Table 2: The results for various quantities for the L\ of fit 10 and sets A, B and C. These 
are described in the text. 
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notation C[ 2 + 2C[ 3 [c5] used in Table E] means that the quantity C[ 2 + 2C[ 3 is determined 
from dg using the L\ as determined from a fit to the same experimental input as used in 
fit 10 in Ref. (2HJ but with the values of L\ and Lg as indicated. The difference between 
the first two determinations of C[ 2 + 2C[ 3 is due to the difference from Fg(Q) from 2Bq. 

The three other sets - A,B, and C - of L\ were generated by different sets of L\ 
and Lg in the agreement region of the two previous sections, Lg « L\ — 0.00035 and 
0.0003 < Z/4 < 0.0006. The extracted constants are generally similar over the three sets 
but still show a dependence on L\ and Lg. There are particular disagreements on values 
obtained via 7^ and 7^. For the former, we can, as we did for the light quark scalar 
radius, attribute the variations over the extracted values to the construction of 7^ on the 
canonical MO solutions. For 7^ s , the argument is similar, as now the relative strength of 
the MO solutions is compensated by the values of Fg s (0)/B and Fg s (0)/B . 

A final conclusion on the value of C[ 2 which is needed for the measurement of V us [26\ 
is not very easy since it depends so strongly on the values of L\ and Lg as well as on which 
dispersive quantities are used to evaluate it. 

5.6 Masses and Decay Constants 

For completeness we present here results for the masses and decay constants for the sets 
of input parameters we use. These can also be found in Tabled We present results for 



Notice that in the corresponding numbers were quoted for the "Main Fit" which used 
the old K e4: data. 

A glance at Table 121 shows that the pion decay constant in the chiral limit Fq can be 
substantially different from the value of about 87 MeV for the case of fit 10. 

We would like to remind the reader that from within SU(3) ChPT all the solutions 
shown in the table are good fits to all fitted experiments. The new sets are preferred under 
the assumptions that the C\ contributing to the scalar form factors and masses vanish as 
well as the constraints from assuming the corrections to the scalar form factors at zero to 
be small. 

6 Conclusions 

In this paper we have computed the full order p 6 expressions in Chiral Perturbation Theory 
in the isospin limit for all pion and kaon scalar form factors. This is described in the first 
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major part of this paper. In addition we have shown various numerical results that allow 
the reader to get a feeling for the size of the various contributions. 

We also found a new relation between the scalar form factors, analogous to the Sirlin 
relation for vector form factors. 

This calculation is important for future determination of the p 6 LEC's and in the 
absence of a better knowledge of these constants the precision on numerical results is 
limited. 

In the second major part of this paper we have used the experimental input used in the 
previous fits of the order p 4 LECs and refitted the L\ as a function of the input values of L\ 
and Lg showing that a rather wide range of these leads to good fits to the included data. 
We also used the solutions of the Muskhelishvili-Omnes problem for the scalar form factors 
with as input the known tttt and KK scattering phases. This allowed to study the scalar 
radii and curvatures as a function of L\ and Lg and obtain consistency relations between 
them. Under the assumptions of vanishing C\ contributions to the scalar form factors at 
t = and the scalar radii these consistency requirements together with the assumption 
that the corrections to the scalar form factors at zero should be small, there is a preferred 
region of L\ and Lg around L\ w 0.00045 and Lg « 0.00015. In Tabled we have shown 
the values of several quantities in this region. 

We have used our results for the curvature to obtain constraints on two order p 6 con- 
stants, C[ 2 and C[ 3 . These are of the expected naive order of magnitude but the precision 
obtained is rather low. 
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A The proof of the symmetry relation 

The relation for the scalar form factors, Eq. (j3.18j) can be proven using the same method 
as used in the second paper of Ref. [27] when the different behaviour of scalar and vector 
currents under charge conjugation is taken in to account. 
We write the Hamiltonian in the form 

H = H + X 3 v 3 . (A.l) 

H is the Hamiltonian in the limit m s = m u . We now do oldfashioned perturbation theory 
in A3, and look at the variation of 

F^'it) (A.2) 
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w.r.t. A3 at A3 = 0. This can be rewritten as an expectation value involving the K° and 
7r~ field, v 3 and su. Applying the operator Gy = Cexp(inV2) with V2 the second generator 
of y-spin and C charge conjugation we obtain that 







<9A, su 



(t) 







A 3 =0 



9A, su 



(t) 



(A.3) 



A 3 =0 



while charge conjugation and the fact that it is a scalar current require the full form factors 
in (jA.3|) to be equal. The first variation w.r.t. A3 of F^ n (t) thus vanishes. 
The same argument can be applied to the combination 



(A.4) 



and in the V^-spin limit both combinations are related. This relation can be brought into 
the form of Eq. (|3.18|) using isospin. 



B Analytical expressions for F$(t) 

For brevity we only write the pion scalar form factor with the light quark densities. The 
expressions for the others can be obtained from the authors. The integrals can be found in 
several places. The functions A, B,. . . can be found in [21] and [21]. The H F ,. . . are defined 
in [21] and the Vi can be found in [3T]. The method to evaluate the Hf was developed in 
[2*T] and the Vi are evaluated as described in [31] using the methods of [41J. 

The result given is not the shortest possible analytical one. There are various relations 
between the integrals that we have not implemented. The reason is that these relations 
involve inverse powers of t and we have used simplifying relations to rewrite all the masses 
in the numerators in terms of m*. and m 2 K . This has as a consequence that, if we had used 
those relations, the numerators of the 1/t poles would not cancel exactly numerically and 
produce possible instabilities for small t. 

We write the pion scalar form factor as 

m*) = 2^0 (V 2) + j- 2 Fs W (t) + jzF%®{t)) . (B.5) 

The superscript indicates the chiral order. The lowest order is simply 

Fg {2) = 1 . (B.6) 

The next order has been calculated in and we agree with their result, it is 

i^ (4) (t) = -lQLlm 2 K - 2ALlml + SL r 4 t - VoL\ml + AL r 5 t + 32L r 6 m 2 K 

+48L r 6 ml + 32L;ml + l/6Z(mJ) - l/2Z(m^) +B(m 2 l ,m 2 v ,t)(l/18ml) 
+B(m 2 K , m 2 K , t)(l/4t) + B(ml, ml, t)(-l/2mj + t) . (B.7) 
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The next order we split in several different parts 

Ff\t) = FUt) + FUt) + FSy(t) + F s \(*) + FSc(t) . (B.8) 

The result for Fg C (t) can also be found in the main text we simply repeat it here for 
completeness. 

F£ c (t) = -128m>^C[ 3 -64m^C[ 5 + 64m>^C[ 6 + 6Amlm 2 K C r 20 

+576mlm 2 K C r 21 + 128m>^C 3 r 2 + 32m 2 CT 2 * + 64m 2 CT 3 * + 16m 2 CT 4 * 
+24m 2 CT 5 * + 32m 2 CT 6 * + 16to 2 C[ 7 * + 16mlC r u t + lOrn 2 ^* - 96m^[ 2 
-128m^C[ 3 - 48m^[ 4 - 64m^C[ 5 - H2m^C[ 6 - 48m^C[ 7 + 144m*C[ g 
+240mlC r 20 + 2A$m*C r 21 + 96mlC r 31 + 128m^C 3 r 2 + 16m 2 K C r 15 t - 6Am 4 K C r w 
+64m^C 2 r + \%2m A K C r 21 - 8C r 12 t 2 - 16C[ 3 * 2 . (B.9) 



Fsn(t) = H F \ml,mlml,ml)(5/6ml) + H F \ml,m 2 K ,ml,ml)(-5/8ml) 
+H F '(ml m 2 , m 2 , m 2 ) (l/18m£) + m 2 , m 2 ,, m 2 ) (m>^) 

+H F '(m 2 K , m 2 K , m 2 , ml ) (-5/ 6m^) 
+H F '(m 2 , m 2 K , m\, ml)(l/2mlm 2 K + 7/8m*) 

i ttF/ / 2 2 2 2\/ 4\ , nF// 2 2 2 2 \ / 4 \ 

+#! {m„,m K ,m K ,m„)(m„) + H 1 (m v ,m K ,m K ,mJ(-mJ 

i rrF// 2 2 2 2\/o 4\ , u-F// 2 2 2 2\/ o /n 4\ 

+H 21 (m ff , m ff , m ff , mj(3mj + if 21 (m ff , m x , m x , toJ(-3/8toJ 
+# 2 F 1 / (m 2 „ m 2 , m^, m 2 ) (3mJ) + ^'(m 2 , m^, m 2 ,, m 2 ) (9/8mJ) 
+H F (ml, ml, ml, m 2 )(2/3m 2 - 2/3*) 

+H F (m 2 n , m\, m\, TO 2 )(-47/24m 2 + 1/6*) + H F (m 2 n , m% m% m 2 )(l/9m 2 ) 
+H F (m 2 K ,ml,m 2 K ,ml){5/Qml + m 2 K - 1/2*) 
+# F (m 2 r ,m 2 r ,m 2 ,m 2 )(-13/9m 2 - 1/8*) 
+# F (m 2 ,,?74,m 2 r ,m 2 r )(119/72m 2 r + l/2m|- + 1/24*) 
+H F (ml,m 2 K ,m 2 K ,ml)(2ml + 1/24*) 
+H[ {m\,m 2 K ,m 2 K ,ml){-2ml - 1/24*) 

+if£ (ml, ml, ml, ml) (6m 2 ) + H F (m 2 , m\, m\, ml ) {-3 /Ami ) 
+H F (m 2 K ,ml,m 2 K ,ml)(Qml) + H F (m 2 ,m 2 K ,m 2 K ,ml)(9/Aml) . (B.10) 



*sv(t) = +V (ml,ml,ml,ml,ml,t,ml)(-mlt + 7/2ml + 2/3t 2 ) 

+V (ml, ml, m\, m\, ml, t, m^)(-ll/12m 2 * + 3/2m* + 1/12* 2 ) 
+V (ml , ml, m 2 , m 2 , m 2 n , t, m 2 )(l/18m^) 
+V (ml,m 2 ,m 2 K ,m 2 K ,ml,t,ml)(l/Amlt - l/&m*) 
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+V (m 2 K , m 2 K , ml, m\, m\, t, m 2 ) (-17/12™^ + 3/2m 4 + l/2t 2 ) 
+V (m 2 K , ml, m 2 K , m 2 , ml, t, m^)(-2/3m^ + 2/3m 4 w + l/8t 2 ) 
+Vo(m 2 rj ,ml,ml,ml,ml,t,ml)(2/9mlml + 13/36m 2 t 

- 7/18m 4 w - l/9m 2 K t - l/12t 2 ) 
+V (m 2 , m 2 , ml, m 2 , m 2 w , t, m^)(l/27m*) 

+V (m 2 , m 2 , m^, m^, m 2 , t, m^)(-5/36m^ + l/ISm* + l/12t 2 ) 
+Vii(m 2 ,m 2 ,m 2 ,m 2 ,m 2 ,t,m 2 J(13/3m 2 t - 4m^ - 4/3t 2 ) 
+\/ n (m 2 , m 2 , m^, , m 2 , t, m 2 n ) (1 1/ Qm 2 J - 3m 4 w - 1/ 6t 2 ) 
+Vu(ml, m 2 , m 2 K , m 2 K , m 2 , t, m 2 )(-l/2m 2 t + 2/3m 4 ) 
+V n (m 2 K , m 2 K , ml, m 2 K , ml, t, m 2 )(53/12m 2 t - 4m 4 n - t 2 ) 
+V 11 (ml,m 2 K ,m 2 K ,m 2 ,ml,t,m 2 r )(5/3mlt - 2ml ~ l/4t 2 ) 
+V 11 (m 2 ,ml,m 2 K ,m 2 K ,ml,t,ml)(-2/9mlm 2 K - 7/9m 2 t 

+ 8/9ml + l/9m 2 K t + l/6t 2 ) 
+V 11 (m 2 ,m 2 ,m 2 K ,m 2 K ,ml,t,ml)(5/9mlt - 1 /3m* - l/6t 2 ) 

iT//2 2 2 2 2 j. 2 \ I ,2\ 

+V 12 {m 7T ,m 7r ,m 7r ,m 7r ,m 7r ,t,m 7r ){-t ) 

+V 12 {m 2 n , ml, m 2 K , m 2 K , m 2 n , t, m 2 )(l/3m 2 t + l/8t 2 ) 

+V 12 (ml,m 2 ,ml,m 2 K ,ml,t,m 2 7r )(-l/6mlt - l/8t 2 ) 

+V 12 (m 2 K ,m 2 K ,ml,m 2 K ,ml,t,ml)(5/Qmlt - l/2t 2 ) 

+V 12 (m 2 K , m 2 K , m 2 K , m 2 , m 2 w , t, m 2 )(l/3m 2 t) 

+V 12 (m 2 , ml, m 2 K , m 2 K , m 2 n , t, m 2 )(-5/18m 2 t + l/9m 2 K t + l/8t 2 ) 

+V 12 {m 2 v ,m 2 v ,m 2 K ,m 2 K ,ml,t,ml){l/9mlt - l/8t 2 ) 

+V 13 {m 2 w , ml, ml, m 2 n , m 2 n , t, m 2 w )(2/3m 2 t - 8ml + 2/3t 2 ) 

+V 13 (ml, ml, m 2 K , m 2 K , ml, t, m 2 )(2m 2 t - 4ml) 

+Vi 3 (ml, m 2 , m 2 K , m 2 K , m 2 , t, m 2 )(-l/2m 2 t + l/3m*) 

+Vi3(ra 2 r ,m 2 r ,m 2 ,m 2 ^,m 2 ,t,m 2 J(13/6m 2 t - 3ml - 5/24t 2 ) 

+Fi 3 (m^, m 2 K , m 2 K , m% ml, t, m 2 )(3/2m 2 t - 2m\. - l/8t 2 ) 

+yiz{rn 2 v ,ml,m 2 K ,m 2 K ,ml,t,ml)(-A/9mlm 2 K - 13/18m 2 t + 7/9m^ 

+ 2/9m 2 K t + l/6t 2 ) 
+l / i 4 (m 2 , ml, ml, ml, ml, t, m 2 )(2/3t 2 ) 
+V M (ml, ml, m 2 K , m 2 K , m 2 n , t, m 2 )(l/2t 2 ) 
+V 1 4(ml,m 2 1 ,m 2 K ,m 2 K ,ml,t,ml)(l/6mlt - l/At 2 ) 
+V u (m 2 K , m 2 K , ml, m 2 K , m 2 n , t, m 2 )(7/24t 2 ) 
+V u (m 2 K , m 2 K , m 2 K , m 2 , m 2 w , t, m 2 )(l/4t 2 ) 

+V u (m 2 , ml, m 2 K , m 2 K , m 2 n , t, m 2 )(-7/18m 2 t + 2/9m 2 K t + l/6t 2 ) 
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ml, ml, ml, ml, ml, t, ml) (4ml ~ 2 / 3t ) 
ml,ml,m 2 K ,m 2 K ,ml,t,ml)(3/2ml - 7/12t) 
ml,m 2 l ,m 2 K ,m 2 K ,ml,t,ml)(-l/2ml + 1/4*) 

2 2 2 2 2 j. 2\/,i 2 j.\ 

m K ,m K ,m 7T ,m K ,m 7T ,t,m n ){Am 7r - t) 

m 2 K ,m 2 K ,m 2 K ,m 2 ,ml,t,ml)(3/2ml - 1/2*) 

m 2 ,ml,m 2 K ,m 2 K ,ml,t,ml)(-l/2ml + 1/12*) 

m 2 ,m 2 ,m 2 K ,m 2 K ,ml,t,ml)(l/2ml - 1/12*) 

mj, m*, m 2 , m 2 , m 2 , t> m ? j (_ 10 / 3m ^ + 4m* + 2/3* 2 ) 

m 2 , m 2 , m 2 ,, m 2 ,, m 2 , *, m 2 )(-ll/12m 2 * + 3/2m* + 1/12* 2 ) 

m 2 , m 2 , m^, m^, m 2 , t, m 2 )(l/4m 2 * - l/2m*) 

m\, m\, ml, m\, m\, t, m 2 )(-3m 2 * + 4m^ + l/2t 2 ) 

m\, m\, m\, m 2 , m 2 n , t, ml)(-mlt + 3/2m^ + l/8t 2 ) 

m 2 v , ml, m\, m\, m\, t, m 2 )(5/12m 2 * - l/2m* - \j\2t 2 ) 

m\, m% m\, m\, m\, t, m 2 w )(-5/12m 2 * + l/2r< + 1/12* 2 ) 

2 2 2 2 2 j 2\/i /oj-2\ 

m 2 , m 2 , m^, m^, m 2 , *, m 2 )(-5/24* 2 ) 

ml,m 2 ,m 2 K ,m 2 K ,ml, t, m 2 ) (1/ 8* 2 ) 

m 2 K ,m 2 K ,m 2 K ,m 2 i ,ml, t, m 2 n ) (- l/8t 2 ) 

m 2 , m 2 , m^, m^, m 2 , *, m 2 )(-l/24* 2 ) 

m 2 , m 2 , m\, m\, m\, t, m 2 )(l/24* 2 ) 

m 2 w , ml, ml, m 2 n , m 2 n , t, m 2 w )(-8/3m 2 * + * 2 ) 

ml,ml,m 2 K ,m 2 K ,ml,t,ml)(~l/3mlt - l/8t 2 ) 

2 2 2 2 2j. 2\/i /oj.2\ 

m 7r ,m ri ,m K ,m K ,m 7r ,t,m 7T )(l/8t ) 

m^, m^, m 2 , m\, m\, t, m 2 )(-2m 2 * + 1/2* 2 ) 

2 2 2 2 2 i 2\t 1 /o 2i\ 

m K , m K , m K , m v , m w , t, mj(-l/2mj) 
m 2 , ml, m 2 K , m\, ml, t, m 2 )(l/3m 2 * - l/8t 2 ) 
m 2 , m 2 , m 2 K , m 2 K , ml, t, m 2 w )(-l/3m 2 * + l/8t 2 ) 
ml, ml, ml, ml, ml, t, m 2 w )(-4/3t) 
m 2 ,m 2 ,m 2 ^,m 2 r ,m 2 ,*,m 2 )(4m 2 -2*) 

2 2 2 2 2 j. 2\/ 2,1 /oj.\ 

m 7r ,m ri ,m K ,m K ,m 7T ,t,m 7T )(-m n + 1/2*) 
m 2 K ,m 2 K ,ml,m 2 K ,ml,t,ml)(4ml - 19/12*) 

2 _2 2 2 2 . ™2\/o 2 



m^,m^,m x ,m^,m 7r ,*,m u )(3m 7r - 5/4*) 
m 2 v ,ml,m 2 K ,m 2 K ,ml,t,ml)(-ml + 1/6*) 
m 2 , m 2 , m 2 , m 2 , m 2 , *, m 2 )(4/3m 2 * - 2/3* 2 ) 
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+V 26 (ml, ml, m 2 K , m 2 K , m 2 w , t, m^)(-2m^ + 4m*) 

+V 26 (ml, m 2 , m\, m\, ml, t, ml)(l/2m 2 n t - m*) 

+V 26 (m 2 K , m^, m 2 , m^, m^, *, m 2 )(-29/12m 2 * + 4m* + 5/24* 2 ) 

+\/ 26 (m^, m\, m\, m\, ml, t, m 2 )(-7/4m 2 * + 3m* + l/8t 2 ) 

+V 2 6(m 2 ,ml,m 2 K ,m 2 K ,ml,t,ml)(5/6mlt - m% - l/6t 2 ) 

+V 27 (ml, ml, ml, ml, ml, t, ml )(-A/3t 2 ) 

+V 27 (ml, ml, m\, m\, ml, t, m 2 )(-* 2 ) 

+V 27 (m 2 n , m 2 , m 2 K , m 2 K , m 2 , t, m 2 )(l/4* 2 ) 

+V 27 (m 2 K , m 2 K , ml, m\, ml, t, m 2 )(-7/12* 2 ) 

+V 27 (m 2 K , m\, m\, m 2 , m 2 n , t, m 2 )(-l/2* 2 ) 

+V 27 (m 2 , ml, m\, m\, ml, t, m 2 )(-l/12* 2 ) 

+V 28 (m 2 n , ml, ml, m 2 n , m 2 n , t, m 2 )(8/3m 2 * - 4/3* 2 ) 

+V 28 (ml, ml, m\, m\, ml, t, m 2 )(-l/2* 2 ) 

+V 28 (ml, m 2 , m\, m\, ml, t, m 2 )(l/2m 2 *) 

+V 28 (m 2 K ,m 2 K ,ml,m 2 K ,ml,t,ml)(-5/Qmlt - 1/12* 2 ) 

+V 28 (m 2 K ,m 2 K ,m 2 K ,m 2 l ,ml,t,ml)(-l/2mlt - l/8t 2 ) 

+V 28 (m 2 ri ,ml,m 2 K ,m 2 K ,ml,t,ml)(l/Qmlt - l/12t 2 ) 

+V 29 (ml , ml, ml, ml, ml, t, m 2 ) (-2/3* 2 ) 

+V 29 (ml, ml, m\, m\, ml, t, m 2 )(-l/2* 2 ) 

+V 29 (ml , m 2 , m\, m\, ml, t, m 2 )(-l/2m 2 * + l/4t 2 ) 

+V 29 (m 2 K , m\, ml, m\, ml, t, ml) (-7 / 2 At 2 ) 

+V 29 (m 2 K , m 2 K , m 2 K , m 2 , m 2 , t, m 2 )(-l/4* 2 ) 

+V 29 (m 2 , ml, m\, m\, ml, t, m 2 )(l/2m 2 * - 1/6* 2 ) 

+V 210 (ml, m 2 n , m 2 n , m 2 n , m 2 n , t, m 2 )(8m 2 - At) 

+V 210 (ml,ml,m 2 K ,m 2 K ,ml,t,ml)(Aml - 2t) 

+V 2l0 (m 2 K ,m 2 K ,ml,m 2 K ,ml,t,ml)(3/2ml - 3/4*) 

+V 2W (m 2 K ,m 2 K ,m 2 K ,m 2 ,ml,t,m 2 7r )(3/2ml - 3/4*) 

+V 2U (m 2 n , m 2 n , ml, m 2 n , m 2 n , t, m 2 ) (-4m 2 * + 8m 4 n ) 

+V211 (ml, ml,m 2 K ,m 2 K ,ml, t, ml ) (-2m 2 * + 4m*) 

+V 21l (m 2 K , m 2 K ,ml,m 2 K ,ml, t, m 2 ) (-3/ 4m 2 * + 3/2m* ) 

+V2ii(m 2 r,^,m 2 ^,m5,m 2 ,*,m 2 )(-3/4m 2 * + 3/2m^) 

+V 2l2 (ml, m 2 n , ml, m 2 w , m 2 w , t, m 2 )(-2* 2 ) 

+V 212 (ml, ml, m\, m\, ml, t, m 2 )(-* 2 ) 

+V 212 (m 2 K , m\, ml, m\, ml, t, m 2 )(-3/8* 2 ) 
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+V 212 (m 2 K , m 2 K ,m 2 K ,m 2 ,ml, t, ml) (-3/ 8t 2 ) 
+V 213 (ml, ml, ml, m\, ml, t, m 2 )(-2t 2 ) 
+V 2 iz{ml, ml, m 2 K , m\, ml, t, ml)(-t 2 ) 
+V 213 (m 2 K , m 2 K , ml, m\, ml, t, m 2 )(-3/8t 2 ) 

+V 213 (m 2 K , m\, m\, m 2 , ml, t, ml )(-3/8t 2 ) . (B.ll) 

FUt) = +^- 2 (-n/81mlm 2 K /m 2 A(ml) - 5/108m^(mJ) + 5/6m^(m 2 ) 

- l/4ml~A(m 2 K ) + 7 / 324ml / m 2 v A(m 2 v ) - 5/54m 2 K A(m 2 ) + m 2 K A{ml) 
+ m 2 K A(m 2 K ) + 16 /81m%/ 'm\ \A{m 2 ^)) 

+ (^^j (-55/162mlm 2 K n 2 - 1015 / ' 432mlm 2 K + 23/108m 2 vr 2 t 

+ 23/18m 2 t - 1171/1296m^7r 2 - 8783/1728m^ + 17/216m^7r 2 i 

+ 17/36m 2 K t - ll/36m^7r 2 - 15/16m^) 
+B(ml, ml, t) 2 (-l/2m 2 t + l/Am* + 2/9t 2 ) 
+B{ml, ml, t)B(m 2 K , m\, t) (-l/24m 2 t + 7/72t 2 ) 
+B(ml, ml, t)B(m 2 v , m% t) (l/27m 2 t + l/ISm*) 
+B(ml, ml, t)B l (ml, ml, t){-m 2 n t + 10/9t 2 ) 
+B(ml , ml, t)B l (m 2 K , m\, t) (-l/12m 2 t + 17/72t 2 ) 

+B(ml , ml, i)S 21 (ml, ml, t) (4/9t 2 ) +jB(m 2 , ml, t)B 2l (m\, m\, t) (l/12t 2 ) 
+B(ml , ml, t)B 22 (ml, ml i)(4/9f) + B(m\, ml, t)B 22 (m\, m\, t)(l/12i) 
+fi(m 2 ,m 2 ,t)(-l/12m 2 A(m 2 ) - 15/4m 2 A(m 2 ) - 2m 2 n A{m 2 K ) 

+ l/9L4(m 2 ) + tA(ml) + 3b/72tA(m 2 K )) 
+B(ml,ml,0)(-l/9mlA(m 2 v ) + l/3m 2 A(m 2 )) 
+B(m 2 K ,m 2 K , t) 2 (l/2tt 2 ) 

+B(m 2 K , m\, t)~B{m\, m\, t){-l/9m 2 n m 2 K + l/24m 2 t - \j21m\t + l/72t 2 ) 

+5(m^, m^, t)B x {ml, ml, t)(-l/12m 2 t + l/4t 2 ) 

+5(m^ m*-, t)B x {m 2 K , m\, t)(5/24t 2 ) 

+B(m 2 K , m\, t)£i(m 2 , m 2 , t)(l/12m 2 t + l/36t 2 ) 

+B{m 2 K ,m 2 K ,t)B 21 ( t)(l/9t 2 ) 

+5(rr4, m^, t)B 2l {m 2 K , m\, t)(l/12t 2 ) 

+B(m^, m^, t)B 22 (ml, ml, f)(l/9t) 

+5(m^, m^, t)B 22 (m 2 K , m\, f)(l/12t) 

+fi(m 2 ^,m 2 r ,t)(-l/4m 2 Z(mJ) - 3/4m 2 A(m 2 ) - l/2m 2 Z(m 2 Y ) 

- 53/360L4(m 2 ; ) + 35/72L4(m 2 ) + 17/120L4(m^)) 
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+B(m 2 K ,m 2 K ,0)(l/18m 2 K A(m 2 v )) 
+B(m 2 rj ,m 2 v ,t) 2 (A/81m 2 7r m 2 K - 7/324m*) 
+B(m 2 , m% *)5i(m 2 , ml, t)(l/27m 2 t) 

(m 2 K ,m 2 K ,t)(l/12mlt - l/27m 2 K t + l/24t 2 ) 
+5(m 2 , m 2 , t)B 2l {m 2 K , m 2 K , t)(l/36t 2 ) 
+B(m 2 , m 2 , t)B 22 {m 2 K , m\, t)(l/36t) 

+B(m 2 ,m 2 ri ,t)(l/l08mlA(m 2 ) + 7/36m^(m 2 ) - 2/9mlA(m 2 K ) 
+ 13/360L4(m^)) 

+B\ml,ml,t)—^(-2mlm 2 K + 14/9m 2 £ - 91/18m* + l9/36m 2 K t 



+ 97/72t 2 ) 
,t 

+ 35/72t 2 ) 



+B e (m 2 K , m 2 K , t)-^-{-mlm 2 K + 103/180m 2 t - 2/3m^ + U/45m 2 K t 



+ J B e (m;,m 2 ,t)^(l/6m 2 t - 5/54m 4 w + l/20m^ + l/24t 2 ) 
+B\{ml , ml, t)^(25/18m 2 t + l7/18m 2 K t - 253/72t 2 ) 
+S e 1 (m 2 Y ,m 2 Y ,t)-^(-71/180m 2 t + 17 /45m 2 K t - 101/72t 2 ) 
+ B\(m 2 ri ,m%t)^- 1 (-l/Qmlt - l/10m 2 K t - 5/24t 2 ) 
+B e 21 (m 2 n , ml, t) ^(55/72^) + ~W 21 {m 2 K , m\, t)-L(13/36t 2 ) 
+5 21 (m 2 , m 2 , t)^- 2 (l/8t 2 ) + 5 22 (m 2 , m 2 , *)-L(-l 5 /2m 2 + 95/36t) 
+^ 2 (m 2 m 2 t) I ^(-3m 2 + 10/9*) 

+ B\ 2 {ml,m%t)^- 2 {-l/2ml + 1/4*) + ^(m 2 , m 2 , t) 2 (4/9t 2 ) 

+£?i(m 2 , m 2 , t)B x [m\, m\, t) (7/36t 2 ) 
+5 1 (m 2 ,m 2 ,t)(l/9L4(m 2 ) + 5/9L4(m 2 ) + 17/18tA(m|-)) 
+Bi(m 2 K , m\, t) 2 (l/12t 2 ) + ^(m^, m^, t)B x {m\, m 2 , t)(l/36t 2 ) 
+B 1 (m 2 Y ,m 2 Y ,t)(23/180£4(m 2 ) - 1 / 12tA(ml) + 2 / 15tA(m 2 K )) 
+B 1 (m 2 ,m 2 ,t)(-l/10tA(m 2 K )) - l/2m^ 2 m 2 K (A(m 2 n )) 2 
+23/18m~H (A(ml)) 2 - A/3m 2 w m K 2 (A(m 2 K )) 2 
+ ll/162m 2 m^/mj (A(m 2 )) 2 - 35/144m 2 /m 2 ; (A(m 2 )) 2 
-7/648m 4 Jm 4 v (A(m 2 )) 2 + 17/36m7^ (A(m 2 K )) 2 
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-7/36m 2 K /m 2 v (A(ml)) 2 - 8/81m 4 K /m 4 (^K))' 

+13/36A(mJ)A(m2) - 7 '/ ' l%A{m 2 v )A{m 2 K ) + 1/36 (A(m 2 v )) 2 -A(m 2 n )A(m 2 K ) 
-1037/144 (A(ml)) 2 - 7/4 (A(m 2 K )) 2 . (B.12) 

FsL(t) = +^(-16/9m>^-16/27m>^3 + 12m>^-24m>^ 

+ 4m^ + 74/9m 4 n L r 2 + 56/27m^ + 16/3m 4 n L r 4 + 16/3m*L r 5 

- 32/3m*L r 6 - 32/3m*L£ + 10A/9m 4 K L r 2 + 86/27m 4 K L r 3 + 8/3m 4 K L r 4 
+ A/3m 4 K L r 5 - 16/3m 4 K L r 6 - 8/3m 4 K L r 8 ) 

+B(ml,ml,t)(8mlm 2 K L r 4 - 16m 2 m 2 ^ - 24m 2 L^t - 8m 2 - 12m 2 

- 12m 2 L^ - 26/3m 2 L^ + 32m 2 L^ + 64/3m 2 L^ + 4&m*L\ + 16m 4 n L r 2 
+ 24m 4 n L r 3 - 36m 4 7r L r 4 - 16m 4 L r 5 + 56m 4 w L r 6 + 24m 4 L r 8 - 32/3m 2 K L r 4 t 

+ 64/3m 2 K L r 6 t) 

+B(ml, ml, 0)(32/3mlm 2 K L r 4 - Q4/3mlm 2 K L r 6 + 16/3m*Ll + 16/3m*L£ 

- 32/3m 4 r L r 6 - 32/3m 4 n L r 8 ) 

+B{m 2 K , m 2 K , t){32mlm 2 K L\ + 8mlm 2 K V 3 - 32m 2 n m 2 K L r 4 - 8m 2 n m 2 K Ll 
+ 32mlm 2 K L r 6 + lQmlm 2 K L r 8 - 4/3m 2 L^ + 8/3m 2 L^ - \6m\L\t 

- Am 2 K L r 3 t - 8/3m 2 K Ut + 16m 2 K L r 6 t + 16/3m 2 K L r s t) 

+B(m 2 K , m 2 K , 0)(4/3mlm 2 K L r A - 8 / 3mlm 2 K L r 6 + 8/3m 4 K L r 4 + 4/3m^L£ 

- lQ/3m 4 K L r 6 - 8/3m 4 K L r 8 ) 
+B(m 2 ,m 2 v ,t)(6A/9mlm 2 K L r 1 + 32/27m 2 m^ - 104/9mlm 2 K L r 4 

- 64/27m 2 m^L£ + 16mlm 2 K L r 6 - 128/9mlm 2 K U 7 + 8/9m 2 L^ 
+ 4/27m 2 n L r 3 t - 4/9m 2 L^ - 2/9m 2 L^ - 16/9m^ 

- 8/27m 4 U 3 + 29/9m 4 L r 4 + 16/27m*L r 5 - 8/3m 4 w L r 6 + 128/9m^ 
+ 40/9m*L£ - Z2/$m 2 K L\t - 16/27m^L^ + 16/9m^t) 

+5 1 (m 2 , ml, t)(-48mlL\t - 16mlL r 2 t - 24m 2 - 8/3m 2 w L r 5 t + 32m 2 Z/ 6 t 

+ Q4/3mlL r 8 t - 32/3m 2 K U 4 t + U/3m 2 K U 6 t + 32L^ 2 + 24L^t 2 

+ \QL r 3 t 2 + 32/3L^ 2 + 16/3L^ 2 ) 
YB x {m 2 K , m\, t)(-32mlL\t - 8m 2 L^t + 44/3m 2 L^ + 4m 2 L^ 

+ 8/3mlL r 6 t - 8m 2 K L\t - 8/3m 2 K L r 5 t + 16m 2 K L r 6 t + 16/3m^L^ + 16L^t 2 

+ 8L^ 2 + 6L r 3 t 2 + 16/3L r 4 t 2 + 4/30t 2 ) 
+B 1 (m 2 l , m 2 , t)(-16/3m 2 L^ - 8/9m 2 L^ + 16/3m 2 L^ + 8/9m 2 L^ 

+ 8/3L^ 2 + 4/3L^t 2 + 8/9L^ 2 ) 
+B 21 (ml,ml,t)(-8L r 1 t 2 - 16L r 2 t 2 - 4U 3 t 2 + 32/3L r 4 t 2 + l6/3L r 5 t 2 ) 
+B 21 (m 2 K , m\, t)(-8L r 2 t 2 - 2L r 3 t 2 + 16/3L^ 2 + 4/3L^t 2 ) 
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+B 21 {m%m%t){-A/3U 2 t 2 - 4/9L r 3 t 2 ) 

+B 22 (ml, ml, t)(327r£L; + 64m 2 Li; + Van^L\ - 16LJt - 32L^ - 8L r 3 t 

+ 32/3L^ + 16/3L£t) 
+B 22 (m 2 K , m 2 K , t)(32mlL r 2 + 8mlL r 3 - 16L r 2 t - AL r 3 t + l6/3L r 4 t + 4/31$) 
+B 22 (m;, m 2 , t)(16/3m 2 L^ + 16/9to=LS - 8/32# - 8/9L£t) 
-128m^m|L^L^ + 1536™^^ + 512mlm 2 K L r 4 L r 8 - 256m 2 m 2 ^) 2 
+256mlm 2 K L r 5 L r 6 - l024mlm 2 K L r 6 L r 8 - 2048m 2 n m 2 K {L r 6 ) 2 + KMm^I^) 
+Z2mlU 1 A{m 2 K ) - 4/3m 2 L£4(m 2 ) + 48m 2 L^4(m 2 ) - 4/9m 2 Lp(m 2 ) 
+52mlLl~A(ml) + 8mlL r 3 A(m 2 K ) - 8/3mlL r 4 A(m 2 v ) - 424/3m^3(m^) 
-172/3mlLlA(m 2 K ) - 4/3m 2 Lp(m 2 ) - 244/3m^^4(m^) - 2^L\A{m\) 
+800/3mlL r 6 A(ml) + 248/ 3ml L r 6 A{m 2 K ) + 32m^Z(mJ) + \§mlU 8 ~A{m\) 
+464/3m 2 L r 8 Z(m 2 ) + A8mlL r 8 A(m 2 K ) - \28mlU 4 L\ + hVlm%L\L\ 
+384r<L^ - 64m^) 2 + 512^1,^ + 384m^ r 5 L r 8 - 64mt(L^) 2 
-1280m* L\L\ - l§8m\{L\) 2 - 512m*(L£) 2 + U/3m 2 K U l A{m 2 ri ) 
+6Am 2 K L r 1 A(m 2 K ) + lQ/3m 2 K L r 2 A(m 2 v ) + 16m 2 K L r 2 A(m 2 K ) + lQ/3m 2 K L r 3 A(m 2 v ) 
+20m 2 K L r 3 A(m 2 K ) - 88/3m 2 K L r 4 A(m 2 ri ) - 104/ '3m 2 K L r 4 A(ml) 
-248 /3m 2 K L\ ~A(m 2 K ) - 32/9m 2 r L^Z(m 2 ) - 52/3m 2 Y L^(m 2 Y ) 
+112/3m 2 ^Z(m 2 ) + 208/3m 2 ^A(m 2 ) + 30A/3m 2 K Ll~A{m 2 K ) 
-Q4/3m 2 K L r 7 A(m 2 v ) + WA/3m 2 K L r 8 A(m 2 K ) - 128m 4 K L r 4 L r 5 + \Q2Am A K L r 4 L r & 
+25Qm%L r 4 L r 8 - 25Qm%(L r 4 ) 2 + 256m%L r 5 L r 6 - 512m%L r 6 L r 8 - 1024m^(L^) 2 
-8/3L[M(m 2 ) - 2AL[tA(ml) - 16L[tA{m 2 K ) - 8L r 2 iA{ml) - 4/9L r 3 tA(m 2 ) 
-12L^Z(m 2 ) - 4L r 3 t~A(m 2 K ) + 4/3L^l(m 2 ) + 20/3U 4 tA(ml) 
+ 16/3L r 4 tA(m 2 K ) + 2/3L^Z(m 2 ) + 10/3L££4(m 2 ) + 10/ 3 L r 5 tA(m 2 K ) . (B.13) 
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